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Abstract. An overview of recent results pertaining to the hydrodynamic description (both New- 
tonian and non-Newtonian) of granular gases described by the Boltzmann equation for inelastic 
Maxwell models is presented. The use of this mathematical model allows us to get exact results 
for different problems. First, the Navier-Stokes constitutive equations with explicit expressions 
for the corresponding transport coefficients are derived by applying the Chapman-Enskog method 
to inelastic gases. Second, the non-Newtonian rheological properties in the uniform shear flow 
(USE) are obtained in the steady state as well as in the transient unsteady regime. Next, an exact 
solution for a special class of Couette flows characterized by a uniform heat flux is worked out. 
This solution shares the same rheological properties as the USE and, additionally, two generalized 
transport coefficients associated with the heat flux vector can be identified. Finally, the problem 
of small spatial perturbations of the USE is analyzed with a Chapman-Enskog-like method and 
generalized (tensorial) transport coefficients are obtained. 
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1. Introduction 

A simple and realistic physical model of a granular system under conditions of rapid flow consists 
of a fluid made of inelastic hard spheres (IHS). In the simplest version, the spheres are assumed to 
be smooth (i.e., frictionless) and the inelasticity in collisions is accounted for by a constant coeffi- 
cient of normal restitution a < 1 [12T]| . At a kinetic theory level, all the relevant information about 
the dynamical properties of the fluid is embedded in the one-particle velocity distribution function 
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/(r, v; t). In the case of dilute gases, the conventional Boltzmann equation can be extended to 
the IHS model by changing the collision rules to account for the inelastic character of collisions 
ifTSl l54l . On the other hand, even for ordinary gases made of elastic hard spheres (a = 1), the 
mathematical complexity of the Boltzmann collision operator prevents one from obtaining exact 
results. These difficulties increase considerably in the IHS case (a < 1). For instance, the fourth 
cumulant 02 of the velocity distribution in the so-called homogeneous cooling state is not exactly 
known, although good estimates of it have been proposed [|22l[65ir76l[80ll . Moreover, the explicit 
expressions for the Navier-Stokes (NS) transport coefficients are not exactly known, but they have 
been approximately obtained by considering the leading terms in a Sonine polynomial expansion 

mmmmmmmmm. 

As Maxwell already realized in the context of elastic collisions [|64l . scattering laws where the 
collision rate of two particles is independent of their relative velocity represent tractable mathemat- 
ical models. In that case, the change of velocity moments of order k per unit time can be expressed 
in terms of moments of order k' < k, without the explicit knowledge of the one-particle velocity 
distribution function. In the conventional case of ordinary gases of particles colliding elastically. 
Maxwell models correspond to particles interacting via a repulsive potential proportional to the in- 
verse fourth power of distance (in three dimensions). However, in the framework of the Boltzmann 
equation one can introduce Maxwell models at the level of the cross section, without any reference 
to a specific interaction potential OTIl . Thanks to the use of Maxwell molecules, it is possible in 
some cases to find non-trivial exact solutions to the Boltzmann equation in far from equilibrium 
situations [l3ll|49l|70l|72l|7l. 

Needless to say, the introduction of inelasticity through a constant coefficient of normal resti- 
tution a < 1, while keeping the independence of the collision rate with the relative velocity, opens 
up new perspectives for exact results, including the elastic case (a = 1) as a special limit. This 
justifies the growing interest in the so-called inelastic Maxwell models (IMM) by physicists and 
mathematicians alike in the past few years [|3lll|5lli|7l[l[IS[IIl[Il[Il[Il[l5l[I3|23l|25l|^ 
|33IMIMl35llMIMl4Ql|4ll|50l[57l|63l6^ Furthermore, it is interesting 

to remark that recent experiments ll56l for magnetic grains with dipolar interactions are well de- 
scribed by IMM. Apart from that, this mathematical model of granular gases allows one to explore 
the influence of inelasticity on the dynamic properties in a clean way, without the need of introduc- 
ing additional, and sometimes uncontrolled, approximations. Most of the studies devoted to IMM 
refer to homogeneous and isotropic states. In particular, the high- velocity tails Il3ll6l l32ll33l[34l[57l 
and the velocity cumulants jSl |33l |42l [62l |63l [68l have been derived. Nevertheless, much less is 
known about the hydrodynamic properties for inhomogeneous situations. 

The aim of this review paper is to offer a brief survey on hydrodynamic properties recently 
derived in the context of IMM and also to derive some new results. Traditionally, hydrodynamics 
is understood as restricted to physical situations where the strengths of the spatial gradients of the 
hydrodynamic fields are small. This corresponds to the familiar NS constitutive equations for the 
momentum and heat fluxes. Notwithstanding this, it must be borne in mind that in granular gases 
there exists an inherent coupling between coUisional dissipation and gradients, especially in steady 
states [|74l . As a consequence, the NS description might fail at finite inelasticity. This does not 
necessarily imply a failure of hydrodynamics in the sense that the state of the system is still charac- 
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terized by the hydrodynamic fields but with constitutive relations more complex than the NS ones 
(non-Newtonian states). We address in this paper both views of hydrodynamics by presenting the 
NS transport coefficients of IMM as well as some specific examples of non-Newtonian behavior. 

The organization of this paper is as follows. The Boltzmann equation for IMM and the mass, 
momentum, and energy balance equations are presented in section[2]Section[3]deals with the first 
few coUisional moments. In sections |4~I - [61 we review the main results referring to the hydrody- 
namic properties of the IMM as a mathematical model of a granular gas. First, in section |4~1 the 
Chapman-Enskog method is applied to the Boltzmann equation and the NS transport coefficients 
are explicitly obtained without any approximation (such as truncation in Sonine polynomial ex- 
pansions). Then, in sections [5] and [6] some shear-flow states where the NS description fails are 
analyzed and their non-Newtonian properties are exactly obtained. In section |T] we consider the 
generalized transport coefficients describing small spatial perturbations about the uniform shear 
flow. Finally, in section [8] the main results are summarized and put in perspective. 



2. Inelastic Maxwell models 

The Boltzmann equation for IMM [[5l[T0l|25l[33l can be obtained from the Boltzmann equation for 
IHS by replacing the term |g • ct| in the collision rate (where g = vi — V2 is the relative velocity of 
the colliding pair and ^ is the unit vector directed along the centres of the two colliding spheres) by 
an average value proportional to the thermal velocity ^j2T/m (where T is the granular temperature 
and m is the mass of a particle). In the absence of external forces, the resulting Boltzmann equation 
is (Si 

(5t + v-V)/(r,v;t) = J[r,v;t|/,/], (2.1) 

where 

J[r,Vi;t|/i,/2] = J^^l^^J da j dv2(«-^ri-l)/i(r,Vi;t)/2(r,V2;t) (2.2) 
is the IMM Boltzmann collision operator. Here, 

n(r,t) = j dv/(r,v;t) (2.3) 

is the number density, is an effective collision frequency, fi^ = 2tx'^^'^ /V{d/2) is the total solid 
angle in d dimensions, and a < 1 is the coefficient of normal restitution. In addition, h is the 
operator transforming pre-coUision velocities into post-collision ones: 

^ 1 + a 

5vi,2 = vi,2 T (g ■ (2.4) 

In Eq. (12.21) the collision rate is assumed to be independent of the relative orientation between the 
unit vectors g and a B31 . In an alternative version [[T0l[IIl[T2l[T3l[T4l[T5l|25l|26l, the collision rate 
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has the same dependence on the scalar product g ■ as in the case of hard spheres. For simplicity, 
henceforth we will consider the version of IMM described by Eqs. (12.11 ) and ( 12.21) . 

The collision frequency u is a free parameter of the model that can be chosen to optimize the 
agreement for a given property between the IMM and the IHS model. In any case, one must have 
oc nT^/"^ to mimic the mean collision frequency of IHS, where 

Tfi f 

defines the granular temperature and V = v — u is the peculiar velocity, 

u(r, t) = J dv v/(r, v; t) (2.6) 

being the flow velocity. 

In a hydrodynamic description of an ordinary gas the state of the system is defined by the 
fields associated with the local densities of mass, momentum, and energy, conventionally chosen 
as n, u, and T. For granular gases, even though kinetic energy is not conserved upon collisions, 
it is adequate to take these quantities as hydrodynamic fields ||29l l30l . The starting point is the 
set of macroscopic balance equations which follow directly from Eq. (12.11) by multiplying it by 
{1, V, and integrating over velocity. These balance equations read 

Dtn + nV ■u = 0, (2.7) 

Dtu + —V ■ P = 0, (2.8) 

mn 

AT + ^ ( V ■ q + P : Vu) = -CT. (2.9) 

an 

In these equations, Dt = dt + u ■ V is the material time derivative, 

P(r,t) = my"dvVV/(r,v;t) (2.10) 



is the pressure tensor. 



is the heat flux, and 



m 

q(r,t) = - / dv\/^V/(r,v;t) (2.11) 



771 f 

is the cooling rate. The energy balance equation (12.91 ) differs from that of an ordinary gas by the 
presence of the sink term —(T measuring the rate of energy dissipation due to collisions. 

The set of balance equations (|2.7| )-( |2.9I) are generally valid, regardless of the details of the 
inelastic model and so their structure is common for both IMM and IHS. It is apparent that Eqs. 
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(|2.7|) - (|2.9I) . while exact, do not constitute a closed set of equations. To close them and get a 
hydrodynamic description, one has to express the momentum and heat fluxes in terms of the hy- 
drodynamic fields. These relations are called constitutive equations. In their more general form, 
the fluxes are expressed as functionals of the hydrodynamic fields, namely 

P = P[n,u,T], q = qKu,T]. (2.13) 

In other words, all the space and time dependence of the pressure tensor and the heat flux occurs 
through a functional dependence on n, u, and T, not necessarily local in space or time. However, 
for sufficiently small spatial gradients, the functional dependence can be assumed to be local in 
time and weakly non-local in space. More specifically, 

Pij =p-r]Ns (^iUj + VjUi - ■ U(5ij^ , (2.14) 

q = -knsVT - lins^n. (2.15) 

In Newton's equation (12.141) and Fourier's equation (12.151) . p = nT = ^Tr P is the hydrostatic pres- 
sure, ?7ns is the shear viscosity, kns is the thermal conductivity coefficient, and /iNs is a transport 
coefficient that vanishes in the elastic case. When the constitutive equations (12.141) and (12.151) are 
introduced into the momentum and energy balance equations, the set (|2.7| )- (|2.9I ) becomes closed 
and one arrives to the familiar NS hydrodynamic equations. 

It must be noticed that, in the case of IHS, the cooling rate C also has to be expressed as a 
functional of the hydrodynamic fields. However, in the case of IMM, ( is just proportional to the 
effective collision frequency u. More specifically ll50l[68l . 



( = (2.16) 

This equation allows one to fix u under the criterion that the cooling rate of IMM be the same as 
that of IHS of diameter a. When the IHS cooling rate is evaluated in the Maxwellian approximation 
[I5l[80l,one gets 

d + 2 

iy=^iyo, (2.17) 

where 

^0= rM^o. ^^'-'J^- (2-18) 

y/7r{d + 2) V m 

The collision frequency is the one associated with the NS shear viscosity of an ordinary gas 
(a = 1) of both Maxwell molecules and hard spheres, i.e., ?7ns = p/i^o = Vo ce = 1. However, 
the specific form (12.181) will not be needed in the remainder of the paper. 

An important problem in monocomponent systems is the self-diffusion process. In that case 
one assumes that some particles are labeled with a tag but otherwise they are mechanically equiv- 
alent to the untagged particles. The balance equation reflecting the conservation of mass for the 
tagged particles is 

+ — V-ji = 0, (2.19) 

nm 
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where xi = rii/nh the mole fraction of the tagged particles and 

]i = mj dVV/i(V) (2.20) 

is the mass flux of the tagged particles. Analogously to the case of Eqs. (I2.7I) - (I2.9I) . one needs a 
constitutive equation for ji to get a closed set of equations. For small spatial gradients, Pick's law 
applies, i.e., 

ji = -/^NsVxi, (2.21) 

where -Dns is the self-diffusion coefficient. 

The Boltzmann equation for IMM, Eqs. (12.11 ) and (12.21) . refers to a monodisperse gas. The 
extension to a multi-component gas is straightforward |l7l [381 Il2l- Instead of a single collision 
frequency u, one has a set of collision frequencies that can be chosen to reproduce the cooling 
rates Qj of IHS evaluated in a two-temperature Maxwellian approximation. In that case, the result 
is El 




u,^ = ^nXi\l^ {— + —], (2.22) 
where aij = (ai + aj)/2 and Tj is the partial granular temperature of species i. 



3. Collisional moments 

As said in the Introduction, the key advantage of the Boltzmann equation for Maxwell models 
(both elastic and inelastic) is that the (collisional) moments of the operator J[f, f] can be exactly 
evaluated in terms of the moments of /, without the explicit knowledge of the latter |[79l . More 
explicitly, the collisional moments of order k are given by a bilinear combination of moments of 
order k' and k — k' with < k' < k. In particular, the second- and third-order collisional moments 
are Ml 

m J dVViV,J[V\f,f] = -uo\2{P^J-pS^,)~CpSiJ, (3.1) 

^ y dV ViV,VkJ[V\f, f] = -z/o|3Q.,fc - [QiSjk + QjSik + qkS^J) • (3.2) 

In Eq. (13. Ih . the cooling rate ( is given by Eq. (12.161) and 

As for Eq. (13.21) . one has 

3^ {l + a)\d-l) 3 
"^1^ = 2^+ 2d{d + 2) -013=2-0,. (3.4) 
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Moreover, 

m 
~2 

is a third-rank tensor, whose trace is the heat flux. In particular, from Eq. (13.21 ) we easily get 



Q,,fc = - \ d^JVy.^kf{Y) (3.5) 



- J dVV^V,J[Y\fJ] = -i/2|ig^. (3.6) 

The evaluation of the fourth-order coUisional moments is more involved and their expressions 
can be found in Ref. Il50l . Here, for the sake of illustration, we only quote the equation related to 
the isotropic moment ifTSllSOll : 

m / dV V^J[Y\f, f] = -1/410M4 + —dY - — %n^i, (3.7) 
J nm nm 

where 

Mi = m j dVyV(V) (3.8) 

is the isotropic fourth-order moment and Iljj = Pij — pSij is the irreversible part of the pressure 
tensor. The coefficients in Eq. (13.71) are 

"^1° = + mr^ ' ^ 

_ il + ay{Ad-l-6a + 3a') 
Ai (3.10) 

Ad{dT^ 

In Eqs. (13.31 ). (13.41 ). and (13.91 ) the collision frequencies z/o|2, i^2\i, and z/4|o have been decomposed 
into a part inherent to the coUisional cooling plus a genuine part associated with the coUisional 
transfers. 

In self-diffusion problems it is important to know the first-order coUisional moment of J[/i, /2] 
with /i 7^ /2. After simple algebra one gets [[38ll42l 

m j dVVJ[V|/i,/2] = -^^o|i(ji-j2), (3.12) 

where 

T^o\i = -1^^ (3.13) 

and js (s = 1, 2) is defined by Eq. (12.201) . 

Before studying the hydrodynamic properties of IMM, it is convenient to briefly analyze the 
so-called homogeneous cooling state (HCS). This is the simplest situation of a granular gas and, 
additionally, it plays the role of the reference state around which to carry out the Chapman-Enskog 
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expansion. The HCS is an isotropic, spatially uniform free cooling state [|2T]| . so the Boltzmann 
equation (12.11) becomes 

dJiv,t) = J[v\fJ], (3.14) 

which must be complemented with a given initial condition /(t>, 0). Since the collisions are in- 
elastic, the granular temperature T(t) monotonically decays in time and so a steady state does not 
exist. In fact, the mass and momentum balance equations (12.71) and (12.81) are trivially satisfied and 
the energy equation (12.91 ) reduces to 

dtT = -CT, (3.15) 



whose solution is given by Haff's law 115511 . namely 

T(0) 



[1 + K(o)^] 



2 ' 



(3.16) 



where we have taken into account that ( oc nT^/"^. The next non-trivial isotropic moment is M^^it). 
The evolution equation for the reduced moment 



Ml{t) 



nm[2T{t)/m] 



(3.17) 



is 



dtM: = - {u,\o - 2C) M; + Ai 



(3.18) 



Note that, while Eqs. (13.151) and (13.161) are valid both for IHS and IMM, Eq. (|3.18|) is restricted to 
IMM. The general solution of Eq. (13.181) is 



M:(t) = [M;(0)-H 



-4(i.4|o/2C-l) 



+ /^4, 



where 



/i4 = 



d^ Ai 



4 1^410 - 2C 

In the one-dimensional case, Eq. (13.91) shows that the difference z/4|o — = —(1 



(3.19) 



(3.20) 



is 

negative definite for any a < 1, so that, according to Eq. (13.191 ). the scaled moment diverges 
with time. On the other hand, if > 2, u^^q — 2(^>0 for any a and, consequently, the moment 
goes asymptotically to the value /i4. In that case, the corresponding fourth cumulant defined by 



a2 = 



4 

M* 

d{d+2r' 



is given by 



6fl 



a] 



0-2 



4rf-7 + 3a(2 



a] 



(3.21) 



(3.22) 
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a 



Figure 1: Plot of the fourth cumulant 02 in the HCS for = 2 (dashed line) and d = 3 (solid line). 



The dependence of a2 on the coefficient of restitution a for d = 2 and 3 is displayed in Fig. [T] 
It is apparent that a2 is always positive and rapidly grows with inelasticity, especially in the two- 
dimensional case. In contrast, the dependence of the IHS 02 on a is non-monotonic and much 
weaker 117611 . 

It has been proven [[T2l[T3l that, provided that /(f , 0) has a finite moment of some order higher 
than two, f{v, t) asymptotically tends toward a self-similar solution of the form 



f{v,t) n 



m 



1 dfi 



2T{t) 



0,(c(t)), c(t) = -=_, (3.23) 

y/2T{t)/m 



where </)/i(c) is an isotropic distribution. This scaled distribution is only exactly known in the 
one-dimensional case tSJ, where it is given by 

All the moments of this Lorentzian form of order higher than two are divergent. This is consistent 
with the divergence of M^(t) found in Eq. ( 13.191) . It is interesting to note that if the initial state 
is anisotropic then the anisotropy does not vanish in the scaled velocity distribution function for 
long times ll50l . As a consequence, while the distribution (|3.24l) represents the asymptotic form 
0/1 (c) for a wide class of isotropic initial conditions, it cannot be reached, strictly speaking, from 
any anisotropic initial state. Whether or not there exists a generalization of (|3.24l) for anisotropic 
states is, to the best of our knowledge, an open problem. 

Although the explicit expression of (ph{c) is not known for d > 2, its high- velocity tail has been 
found to be of the form El [331 [53 

0,,(c)~c-'^-^("), (3.25) 
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where the exponent s{a) is the solution of the transcendental equation 



1 - 



a 



2F1 



s 1 d 3 + 2a - 
2' 2' 2' 4 



l + aV£(3W|) 



(3.26) 



2-Fi[a, h] c; z] being a hypergeometric function |[T1. Equation (13.251) implies that those moments of 
order k > s{a) are divergent. 

The evolution of moments of order equal to or lower than four for anisotropic initial states has 
been analyzed in Ref. Il50l . 



4. Navier-Stokes hydrodynamic description 

The standard Chapman-Enskog method fTf\ can be generalized to inelastic collisions [|2T1l to ob- 
tain the dependence of the NS transport coefficients on the coefficient of restitution from the Boltz- 
mann equation [HI HH |45l gSl ED |5l and from the Enskog equation 031 IS ill [60l|. Here the 
method will be applied to the Boltzmann equation (12.11) for IMM. 

In order to get the hydrodynamic description in the sense of Eq. (12.131) . we need to obtain 
a normal solution to the Boltzmann equation. A normal solution is a special solution where all 
the space and time dependence of the velocity distribution function takes place via a functional 
dependence on the hydrodynamic fields, i.e., 

/ = /[v|n,u,T]. (4.1) 

This functional dependence can be made explicit by the Chapman-Enskog method if the gradients 
are small. In the method, a factor e is assigned to every gradient operator and the distribution 
function is represented as a series in this formal "uniformity" parameter, 

/ = f'^ + 6/« + e^f^'^ + ■■■. (4.2) 

Insertion of this expansion in the definitions of the fluxes (12.101) and (12.1 II) gives the corresponding 
expansion for these quantities. Finally, use of these expansions in the balance equations (I2.7I) - (|2.9I) 
leads to an identification of the time derivatives of the fields as an expansion in the gradients, 

d, = df^ +edP +^d?^ + ■■■ . (4.3) 

The starting point is the zeroth order solution. The macroscopic balance equations to zeroth 
order are 

df^n = Q, af^u = 0, df^T = -CT. (4.4) 

Here, we have taken into account that in the Boltzmann operator (12.21) the effective collision fre- 
quency V oc nT^/'^, and hence the cooling rate C is a functional of / only through the density n and 
granular temperature T [see Eq. (12.161) 1. Consequently, (^^^^ = C,. To zeroth order in the gradients 
the kinetic equation (12.11) reads 

5*/(°) = J[V|/(°),/(°)]. (4.5) 
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This equation coincides with that of the HCS, Eq. (13.141) . Moreover, since f^'^^ must be a normal 
solution, its temporal dependence only occurs through temperature and so it is given by the right- 
hand side of Eq. (13.231 ). except that n — )■ n(r, t) and T T(r, t) are local quantities and v — 
V = V — u(r, t). The normal character of f^^^ allows one to write 



so that Eq. (14.51) becomes 



5(o)^(o)(v) = -CT^/W(V) = 1^ . V/(°)(V), (4.6) 



^^■V/W(V) = J[V|/W,/(°)]. (4.7) 



Since /^^^ is isotropic, it follows that 



ij — y^tj 

Therefore, the macroscopic balance equations to first order give 



P!i^=pS^., q(°^=0. (4.8) 



mn^ d 



D\''n = -nV -u, D);'n = ~-^, Df^T = — ■ u, (4.9) 



where Df^ = d^^^ + u ■ V. To first order in the gradients, Eq. (12.11 ) leads to the following equation 
for/W: 

(Sf + C) /(i)(V) = - + V ■ v) /(°)(V), (4.10) 
where C is the linearized collision operator 

I da I dV2 (a-^ri - l) [/(°)(Vi)/(^nV2) + /(°nV2)/(^HVi)] • 



(4.11) 



Using (14.91) . the right-hand side of Eq. (14.101) can evaluated explicitly, so the integral equation 
for /(^) can be written as 

(^(9f^ + £) f'-^HV) = A(V) ■ VlnT + B(V) ■ Vlnn + C(V) : Vu, (4.12) 

where 

i^^f'"') - -l^f'"'^ (4.13) 

B ^ -V/(°) - (4.14) 
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The structure of Eq. (14.121 ) is identical to that of IHS, except for the detailed form (14.1 II ) of the 
linearized Boltzmann collision operator L. In spite of the advantages of IMM, Eq. (14.121 ) is math- 
ematically rather intricate and its solution is not known. In the case of IHS, a trial function (based 
on a truncated Sonine polynomial expansion) for /*^^^ is proposed. The coefficients in the trial 
function, which are directly related to the NS transport coefficients, are obtained in an approxi- 
mate way by taking velocity moments. On the other hand, the use of a trial function is not needed 
in the case of IMM and the NS transport coefficients can be obtained exactly. The key point is that, 
upon linearization of Eqs. (13.11) and (13.61 ). one has 

m y ^YV,V^Cf'^\\) = uo\2P^p, (4.16) 

^ I dVy2V£/«(V) = z/2|iq(^). (4.17) 
Now we multiply both sides of Eq. (|4.12|) by mViVj and integrate over V to obtain 

(af + u^\2)P^^ = -P (v,u, + V,n, - ■ . (4.18) 

This equation shows that Pj^p is proportional to the right-hand side divided by a collision fre- 
quency. Therefore, P.l^p oc p/nT^/'^ = T^/^ and so 

5(o)p(i) = _^p(i), (4.19) 

where we have taken into account Eq. (14.41) . As a consequence, the solution to Eq. (14.181) is 

P-P = -ms (v,n, + V,n, - ■ u6,,^ , (4.20) 

where 

Vns = (4.21) 

Comparison with Eq. ( 12.141) allows one to identify Eq. ( 14.211 ) with the NS shear viscosity of IMM. 

Let us consider next the heat flux. Multiplying both sides of Eq. (14.121 ) by ^mV'^\ and inte- 
grating over V we get 

+ q(^) = -^(1 + 2a,)^Vr- ^a,^Vn, (4.22) 
V / 2 m 2 m 

where use has been made of Eq. (14.171) . Here, 02 is the fourth cumulant of f^^\ whose expression 
is given by Eq. (13.221) . The right-hand side of Eq. (14.221) implies that the heat flux has the structure 

q(i) = -kns VT - ^iNS Vn. (4.23) 



12 



V. Garzo and A. Santos 



Hydrodynamics of inelastic Maxwell models 



By dimensional analysis, kns T^/^ and /iNs oc T^/^. Consequently, 



C 



3 T 
2knsVT + ( -/iNs + '^Ns— ) Vra 



(4.24) 



where in the last step we have taken into account that C oc nT^^"^. Inserting this equation into Eq. 
(I4.22|) . one can identify the transport coefficients as 

pd + 2 l + 2a2 

i^Ns = ^, (4.25) 

m 2 V2\i - 2C 

/iNS = -T— 7^ T7- (4.26) 

n 1 + 2a2 U2\i - f C 

Equations (14.211) . (14.251) . and (14.261) provide the NS transport coefficients of the granular gas 
modeled by the IMM in terms of the cooling rate C,, the collision frequencies uq\2 and U2\i^ and the 
HCS fourth cumulant a2. Making use of their explicit expressions, Eqs. (12.161) . (13.31) . (|3.4I) . and 
(13.221) . respectively, the a-dependence of the transport coefficients is given by 

^^^ = ^°(l + a) [?>d + 2 + {d-2)aY ^^'^^^ 

8(rf-l)[5 + 4rf-9a(2-c.)] 
'^\l + a){d-A + 3da) [4d - 7 + 3a(2 - a)] ' ^ ^ 

T 16(1 - a) [2d^ + 8d - 1 - 6{d + 2)a + 9a^] 
^""^ ~ '"V (1 + a)2 (d - 4 + 3da) [4d - 7 + 3a{2 - a)] ' ^ ^ ^ 

where ?7o = (rf + 2)(p/2z/) = p/uo and = [d{d + 2) / 2{d — l)]{rio / m) are the NS shear viscosity 
and thermal conductivity coefficients in the elastic limit (a = 1), respectively. 

It is interesting to rewrite Eq. ( 14.231) using T and nT^/^ as hydrodynamic variables instead of 
T and n [|5n[66l . In fact, both the collision frequency u and the cooling rate are proportional to 
nT^f'^. In these variables, the heat flux becomes 

q(i) = -<sVT - /iNsT-^/^V [nT^'^) , (4.30) 

where k^s = '^ns — fJ'Nsin/2T). Inserting Eqs. (14.251) and (14.261) . we easily get 



K 



_ p d + 2 1 + |a2 
" m 2 Z.2I1-IC 



8 [1 + 2rf - 3a{2 - a)] 
°(l + a)2[4rf-7 + 3«(2-«)]' ^ ^ 

The presence of the term d — A + 3da in the denominators of Eqs. (14.281) and ( 14.291 ) implies 
that the heat flux transport coefficients kns and /ins diverge when a tends (from above) to ao = | 
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and ao = I for (i = 2 and d = 3, respectively |[T9l[68l- However, the coefficient is finite for 
any a and any d > 1. 

The one-dimensional case deserves some care. As is known, the thermal conductivity in the 
elastic limit, kq, diverges at d = 1 [|67l . Surprisingly enough, the thermal conductivity is well 
defined at d = 1 for inelastic collisions (a < 1). Taking the limit li — )■ 1 in Eq. (14.251) one gets 
kns = (18p/mz/) /(I — a^). On the other hand, the coefficient /x^fs vanishes at a = 1 but diverges 
for a < 1 if d = 1. 




a 



Figure 2: Plot of the reduced NS shear viscosity ?7ns/'7o for d = 2 (dashed line) and d = 3 (solid 
line). 



Let us try to understand the origin of the divergence of kns and /ins when a — )• ao for d = 2 
and 3. The integral equation (14.121) suggests that its solution has the form 

/(i)(V) = A{V)V ■ VT + B{V)V ■ Vn + C{V) (v^, - ^^'5,,) V,u„ (4.32) 

where A{V), B{V), and C{V) are unknown functions that only depend on the magnitude of veloc- 
ity. The solvability conditions of Eq. (14.121 ) imply that 

j dVV^AiV) = J dW^B{V) = 0. (4.33) 

The transport coefficients are directly related to velocity integrals of ^(V^), B{V), and C{V). 
Specifically, 

K^s = -—j dYV^A{V), (4.34) 

Tfi f 

fi^s = -i^ dVV^B{V), (4.35) 
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a 



Figure 3: Plot of the reduced NS heat flux coefficients kns/^o (solid lines), ji-^^ / {kqT / n) (dashed 
lines), and k'^^/kq (dash-dotted lines) for d = 2 and (i = 3. The quantities k^^ and i^i^^ diverge at 
ao = I (c? = 2) and = | ((i = 3). 




a 

Figure 4: Plot of the reduced NS self-diffusion coefficient Z^ns/^o for any dimensionality d. 
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in I 

The corresponding expression for the modified thermal conductivity coefficient k'^^ is analogous 
to Eq. (14.341) . except for the replacement ^(V) -)> A'(y) = ^(V) - S(V)(n/2T). Equations 
(I4.32I) - (I4.36I) are formally valid for both IMM and IHS. In the former case, however, the algebraic 
high-velocity tail /(°) ~ 't^-'^-^C") [cf. Eq. (l3^ 1 implies, according to Eqs. (03]) - (l4rT5]) . that 
A ~ B ~ Y~d~s{a)+i Q ^ y-d-s{a)^ could thcrcforc expect that the unknown functions 
defining f^^^ also present algebraic tails of the form ^ ~ i3 ~ y-d-a(a)^ ^ y-d-a {a) ^ 
C ~ The convergence of r^Ns and of k'^^ implies that c{a) > 4 and a' (a) > 4 for all a 

and d. However, the divergence of kns and /ins at a = leads to a{a) < 4 for a < if d = 2 
or d = 3. This means that, although /'^^''(V) is well defined for any a, its third-order velocity 
moments (such as the heat flux) might diverge due to the high-velocity tail of the distribution. This 
singular behavior is closely tied to the peculiarities of the IMM since the high-velocity tail of /^^^ 
in the case of IHS is exponential [|20l [37l [801 rather than algebraic. 

To close the evaluation of the NS transport coefficients, we now consider the self-diffusion 
coefficient defined by Eq. (12.211) . It is given by ll42l 



n - ^ 



m - K 

= ^^ n 1 (4.37) 

(1 + a)^ 

where Dq = dp /u is, the self-diffusion coefficient in the elastic limit and we have made use of 
Eqs. (12.161) and (13.131) in the last step. In contrast to the other transport coefficients, the reduced 
self-diffusion coefficient D^s/Dq is independent of the dimensionality of the system. 

Figures [2l-l4] depict the a-dependence of the reduced NS transport coefficients ?7ns/'7oj f^Ns/f^o, 
l-i'Ns / {f^oT / n) , k'^^/kq, and D^s/Dq for d = 2 and d = 3. All of them increase with increasing 
dissipation. As for the influence of d, it depends on the transport coefficient under consideration. 
While, at a given value of a, the shear viscosity increases with dimensionality, the opposite happens 
for the heat flux coefficients. This is especially apparent in the cases of and ^uns since their 
divergence occurs at a smaller value a = ao for c? = 2 than for = 3. Finally, as said above, 
the reduced self-diffusion coefficient is independent of the dimensionality. It is noteworthy that 
the coefficient ji^^/ {kqT /n), which vanishes in the elastic case, becomes larger than kns/'^o for 
sufficiently high inelasticity. 



5. Uniform shear flow 

The hydrodynamic description in the preceding section applies to arbitrary degree of dissipa- 
tion provided that the hydrodynamic gradients are weak enough to allow for a NS theory. The 
Chapman-Enskog method assumes that the relative changes of the hydrodynamic fields over dis- 
tances on the order of the mean free path are small. In the case of ordinary fluids this can be 
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controlled by the initial or boundary conditions. For granular gases the situation is more com- 
plicated, especially in steady states, since there might be a relationship between dissipation and 
gradients such that both cannot be chosen independently. In spite of the above cautions, the NS ap- 
proximation is appropriate in some important problems, such as spatial perturbations of the HCS 
for an isolated system and the linear stability analysis of this state. In this section and the two 
next ones we obtain in an exact way some non-Newtonian hydrodynamic properties of a sheared 
granular gas modeled by the IMM. 

The simple or uniform shear flow (USF) state is one of the most widely studied states, both for 
ordinary [l49| and granular gases ll24l [53l . It is characterized by a constant density n, a uniform 
granular temperature T, and a linear velocity profile Ux = ay, where a is the constant shear rate. 
At a microscopic level, the USF is characterized by a velocity distribution function that becomes 
uniform in the local Lagrangian frame, i.e., 

/(r,v;t) = /(V,t). (5.1) 

In this frame, the Boltzmann equation ( 12.11 ) reduces to 

d,f{Y) - aVy^JiV) = J[V|/, /]. (5.2) 

Equation (|5.2I) is invariant under the transformations 

{V,,Vy)^{-V,,-Vy), (5.3) 

Vj^-Vj^ Vj^Vk, j,k^x,y. (5.4) 

This implies that if the initial state /(V, 0) is consistent with the symmetry properties (15.31) and 
( 15.41) so is the solution to Eq. (15.21 ) at any time t > 0. Even if one starts from an initial condition 
inconsistent with (15.31) and ( 15.41) . it is expected that the solution asymptotically tends for long times 
to a function compatible with (15.31) and (15.41) . 

The properties of uniform temperature and constant density and shear rate are enforced in 
computer simulations by applying the Lees-Edwards boundary conditions ||49l [59l , regardless 
of the particular interaction model considered. In the case of boundary conditions representing 
realistic plates in relative motion, the corresponding nonequilibrium state is the so-called Couette 
flow, where density, temperature, and shear rate are no longer uniform iTTTll . 

According to the conditions defining the USF, the balance equations (12.71) and (12.81) are satisfied 
identically, while Eq. (12.91) becomes 

dtT = -^Pxya - (T. (5.5) 

This balance equation shows that the temperature changes in time due to two competing effects: the 
viscous heating term —Pxyd and the inelastic coUisional cooling term CT. Depending on the initial 
condition, one of the effects prevails over the other one so that the temperature either increases 
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or decreases in time. Eventually, a steady state is reached for sufficiently long times when both 
effects cancel each other. In this steady state 

-^ = ^i (5.6) 

p 2 a 

This relation illustrates the connection between inelasticity (as measured by the cooling rate 0' 
irreversible fluxes (as measured by the shear stress P^y/p), and hydrodynamic gradients (as mea- 
sured by the shear rate a). 

The rheological properties are related to the non-zero elements of the pressure tensor consistent 
with Eqs. (15.31) and (15.41) . namely Pxy, Pxx, Pyy, and p. The remaining d — 2 diagonal elements are 
equal, by symmetry, so that P,^ = ■ ■ ■ = P^^ = i^dp — Pxx — Pyy) /{d — 2). In order to obtain these 
four independent elements, we complement Eq. (15.51) with the equations obtained by multiplying 
both sides of Eq. (15.21) by {V^Vy, V^, Vy} and integrating over velocity. The result is 

dtPxy + O-Pyy = —^0\2Pxyj (5.7) 

dtPxx + 2aPxy = -uo\2 (Pxx -p)- Cp, (5.8) 

dtPyy = -J^0\2 {Pyy - P) - Cp, (5.9) 

where we have taken into account Eq. (13.11 ). 
5.1. Steady-state solution 



The steady-state solution of Eq. (15.91) is simply 

c 

p* =1 — ^ 

d 1 + a 
~ 2d + l- a' 



(5.10) 



Here we have introduced the reduced pressure tensor P-* = Pij/p. Substitution into Eq. (15.71) 
yields, again in the steady state. 



p* 

xy 



Next, Eq. (15.81) gives in the steady state 



a d^{d + 2) 
'u2{d+l-a)^' 



(5.11) 



PL = 1- — + 2^(1-^ 

^^012 Z^o|2 V ^^012 

_ d 1 + a a^ d^{d + 2f 

2d+l-a ^^(d+l-a)3(l + a)' ^" 
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Equations (15.1 II ) and (15.121) are not closed since the ratio u/a, which yields the steady-state tem- 
perature for a given shear rate a, must be determined. This is done by elimination of P*^ between 
Eqs. (1531 ) and ( 15.111 ) with the result 



1^ 2 C 



{I - a^){d + I - af 
Using this result in Eqs. (15.111) and (15.121) . we obtain the a-dependence of P*y and P* 



(5.13) 



d^{d + 2){l-a'^) 
2v^(d + 1 - a) 



^4 = - r^., -■ (5.14) 



^d d + ?>- {d+l)a 
/ d + 1 — a 

Equations (15.101) . (15.141) . and (15.151) provide the explicit expressions of the relevant elements of 
the reduced pressure tensor as functions of the coefficient of restitution a and the dimensionality 
d. Since u"^ oc T, Eq. (15.131) shows that the steady-state temperature is proportional to the square 
of the shear rate. Equation (15.131) can also be interpreted as expressing the reduced shear rate a/u 
as a function of a. As a consequence, no matter how large or small the shear rate a is, its strength 
relative to the stationary collision frequency u is fixed by the value of a, so one cannot choose 
the steady-state value oi a/u independently of a. It is important to notice that, according to Eqs. 
(ISTTOl) and (15131) . P*^ + {d- l)P*y = d. This implies that P*^ = P*y, even though the directions 
y and z are physically different in the geometry of the USE. 

In order to characterize the rheological properties in the USE, it is convenient to introduce a 
generalized shear viscosity r] and a (first) viscometric function ^ by 



Pxy = -V^, (5.16) 



dux 
dy 

2 



Pxx-Pyy = ^(^^j . (5.17) 

The second viscometric function vanishes as a consequence of the property Pzz = Pyy. Erom Eq. 
(I5.11|) one obtains 



where we recall that rjQ = {d + 2){p/2u) = p/uq is the NS shear viscosity at a = 1. Analogously, 
Eqs. (I5J01) and (l5J2l) yield 

2 r d 



^^; = ^, I , (5.19) 

1 + a \d + 1 — a J 

where \E'o = 2?7q/p is the corresponding Burnett coefficient in the elastic limit [|27l . 
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a 



Figure 5: Plot of the reduced rheological functions 77/770 (solid lines) and ^/^^o (dashed lines) in 
the steady-state USF for d = 2 and d = 3. 




0.0 0.2 0.4 0.6 0.8 1.0 
a 

Figure 6: Plot of the reduced shear rate a/i^o in the steady-state USF for d = 2 (dashed line) and 
d = 3 (solid line). 
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a 



Figure 7: Plot of the fourth cumulant a2 in the steady-state USF for d = 3. This quantity diverges 
at etc - 0.046. 

Figure [5] shows the a-dependence of the rheological quantities 77/779 and \l//\l/o for d = 2 and 
d = 3. It is apparent that 7//7/0 is a monotonically decreasing function of inelasticity, this effect 
being more pronounced in the two-dimensional case than in the three-dimensional one. This de- 
crease contrasts dramatically with the behavior of the NS shear viscosity, as seen in Fig. [2l This 
confirms that the transport properties in the steady-state USF are inherently different from those of 
the NS description [17411 . Another non-Newtonian feature is the existence of normal stress differ- 
ences in the shear flow plane. What is interesting is that the viscometric coefficient ^ measuring 
this effect strongly deviates (in general) from its elastic Burnett-order value ^q- We observe from 
Fig.[5]that this effect is again more significant for = 2 than for d = 3. Moreover, in the former 
case ^/^o monotonically decreases with decreasing a, while it reaches a minimum at a = 0.25 
in the three-dimensional case. To complement this discussion, it is worth plotting the steady-state 
reduced shear rate a* = a/u^ versus a. This is done in Fig. [6] for d = 2 and = 3. Since 
a = dux/dy is the only gradient present in the USF, the ratio a/u^ measures the relative strength 
of the hydrodynamic gradients and thus the departure from the homogeneous state. Therefore, it 
plays the role of the Knudsen number. Figure [6] shows that a* increases with inelasticity, having 
an infinite slope at a = 1 . The influence of dimensionality on this quantity is much weaker that in 
the cases of 77/770 and \l//\E'o- 

Although all the previous results in this section are exactly derived from the Boltzmann equa- 
tion (15.21) for EMM, the solution to this equation is not known. However, we can get some indirect 
information about the distribution function /(V) through its moments. In principle, the hierarchy 
of moment equations stemming from Eq. (15.21) can be recursively solved since the equations for 
moments of order k involve only moments of the same and lower order. Equations (15.101 ) and 
(I5.12|) - (|5.15l) give the second-order moments. The next non trivial moments are of fourth-order. 
They were obtained (for d = 3) in Ref. [1731 as the solution of a set of eight linear, inhomogeneous 
equations. The results show that the fourth-order moments are finite for a > olc, where olc — 0.046 
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is a critical value below which the fourth-order moments diverge. This implies that the distribution 
function exhibits a high-energy tail of the form /(V) ~ V~'^~^'^"\ so that s'{a) < 4 if a < ac for 
d = 3. This tail in the USF is reminiscent of that of the HCS [see Eq. (13.251 )1. On the other hand, 
the fourth-order moments are finite in the HCS for > 2 and any value of a [see Eq. (13.221) 1. This 
suggests that s'{a) < s{a), i.e., the shearing enhances the overpopulation of the high-velocity tail. 
As an illustration. Fig. |7] displays the fourth cumulant 02 of the USF, defined by Eq. (|3.21l) . as a 
function of a for d = 3. Comparison with Fig.[T]shows that this quantity is much larger in the USF 
than in the HCS. 

5.2. Unsteady hydrodynamic solution 

The interest of the USF is not restricted to the steady state. In general, starting from an arbitrary 
initial temperature T(0), the temperature T(t) changes in time according to Eq. (15.51) either by 
increasing (if the viscous heating term dominates over the coUisional cooling term) or decreasing 
(in the opposite case). After a short kinetic stage (of the order of a few mean free times) and 
before reaching the steady state, the system follows an unsteady hydrodynamic regime where the 
reduced pressure tensor P*j{t) depends on time through a dependence on the reduced shear rate 
a*{t) = a/z/o(t), in such a way that the functions P*Aa*) are independent of the initial condition 



mm. 



Taking into account that z/q oc T^/^, one has 



dta* 




(5.20) 




(5.21) 



Insertion of this property into Eqs. (I5.7| )- (|5.9I) yields 




(5.22) 




2 



p* 



K\2 {p;y - 1) + c 



(5.23) 



a* 



yy /** _|_ 2a^ p* 

^ d xy 




^o|2(^;.-l) + C + 2a*P;^ 



(5.24) 



* _|_ 2al p* 
d xy 



where we have called v^. 



0|2 



1^012/1^0 and C = C/z/Q. 
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Equations (15.221) and ( 15.231 ) constitute a set of two coupled nonlinear first-order differential 
equations for the elements P*^ and P*y. Their numerical solution, with appropriate boundary con- 
ditions 11741 . provides the hydrodynamic functions P*y{a*) and P*^{a*). Moreover, it is straightfor- 



yy^ 



ward to check that Eqs. (15.231) and (15.241) are consistent with the relationship P*^ + {d — l)Pyy = d. 



This means that the knowledge of Pyy{a*) suffices to determine P*x{a*) and that P*^{a* 



yy^ 



This generalizes the analogous relations (in particular, a vanishing second viscometric 



function) obtained above in the steady state. Once P^ 



xy\ 



mdP;y{a* 



are known, the shear-rate 



dependence of the generalized shear viscosity r] and (first) viscometric function \E', defined by Eqs. 
(15.161) and (15.171) . can be obtained. These two functions are plotted in Fig.[8]for a = 0.6, a = 0.8, 
and a = 1 in the three-dimensional case. The top panel clearly shows that the shear viscosity 
exhibits shear thinning, i.e., it decays with increasing reduced shear rate. As a* increases the in- 
fluence of inelasticity on the shear viscosity becomes less important. The top panel of Fig. [8] also 
shows that the NS value (i.e., the value at a* = 0) of the shear viscosity increases with increasing 
inelasticity, in agreement with Fig.|2l On the other hand, the steady-state values (which correspond 
to different values of a*) decrease as inelasticity increases, in agreement with Fig. |5l Analogous 
features are presented by the viscometric function plotted in the bottom panel. 



Although the determination of P.L{a* 



and Pyy{a* 



involves numerical work, one can obtain 



analytically those functions in the vicinity of the steady state by means of the derivatives da*P*y 
and da-Pyy evaluated at the steady state. This requires some care because the fractions on the right- 
hand side of Eqs. (15.221) and (15.231) become indeterminate in the steady state since the numerators 
and the denominator vanish identically. This difficulty can be solved by means of L'Hopital rule 
[|39l . Therefore, in the steady-state limit Eqs. (15.221 ) and (15.231 ) become 



d ,p* = 

" '^y a* 



P: 



xy 



P* _|_ n*f) ^ P* 
yy ' ^a*J- yy 



d *P 



yy 



xy 

2 

a* 



P*y + a*da'P, 



xy 



p* 



yy 



d ^0\2^°-*^yy 
9 P* _|_ n*f) . P* 



(5.25) 
(5.26) 



Elimination of da*Pyy gives a cubic equation for da*P*y, 



xy 



a*' {da^p:yy + 2du*^,a* {da*p:yY 

+^dP:yP;y + {dp;y - 2p;/) (p;, + d'^^ = o, 



* 2 

d'^-^+Sa* {dP;y-P:y') 



d *P* 

Ua*J- xy 



(5.27) 



with coefficients that are known functions of a. The real root of Eq. (15.271) gives the physical 
solution. From it we simply get 



d .P 



yy 



2P* — 

^^a*P4 + dz/, 



P* + n*d , P* 
^ xy~ ^ ^a*^ xy 



0|2 



+ a*29„.P* 



(5.28) 



xy 
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Figure 8: Plot of the reduced generalized shear viscosity 77/770 (top panel) and of the reduced 
generalized viscometric function ^ /'^o (bottom panel) versus the reduced shear rate a/uQ in the 
unsteady USF for d = 3 and three coefficients of restitution: a = 0.6 (dotted curves), a = 0.8 
(dashed lines), and a = 1 (solid lines). The circles denote the steady-state values. 



6. Couette flow with uniform heat flux. LTu flow 

The planar Couette flow corresponds to a granular gas enclosed between two parallel, infinite 
plates (normal to the y axis) in relative motion along the x direction, and kept at different temper- 
atures. The resulting flow velocity is along the x axis and, from symmetry, it is expected that the 
hydrodynamic fields only vary in the y direction. 

Despite the apparent similarity between the steady planar Couette flow and the USF, the former 
is much more complex than the latter. In contrast to the USF, the temperature is not uniform and 
thus a heat flux vector q coexists with the pressure tensor Pij iTTTll . In general, inelastic cooling 
and viscous heating are unbalanced, their difference dictating the sign of the divergence of the heat 
flux ir77l[83l . More explicitly, the energy balance equation (12.91) in the steady state reads 



dqy d 



= -CnT + P^ya, (6.1) 
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where we have again called a = dux/dy. However, in contrast to the USF, the shear rate a is not 
uniform, i.e., the velocity profile is not linear. Conservation of momentum implies [see Eq. ( 12.81) 1 
Pxy = const and Pyy = const. 

The key difference between the balance equations (15.61) and (16.11) is the presence of the diver- 
gence of the heat flux in the latter. Therefore, Eq. (16.11) reduces to Eq. (15.61) if dyqy = 0, even if 
Qy ^ and dyT ^ 0. This yields a whole new set of steady states where an exact balance between 
the viscous heating term and the coUisional cooling term occurs at all points of the system. This 
class of Couette flows has been observed in computer simulations of IHS and studied theoretically 
by means of Grad's approximate method and a simple kinetic model lf8Tl[82ll . Interestingly, an ex- 
act solution of the Boltzmann equation for IMM supports this class of Couette flows with uniform 
heat flux iTTSl. 

In the geometry of the planar Couette flow, the Boltzmann equation (12.11) becomes 

VydJ = —J[f,f], (6.2) 

where we have particularized to the steady state and have introduced the scaled variable s as 

ds = ujoi2dy, (6.3) 

where ^ 

^0, - -0, - C = (6.4) 

An exact normal solution of Eq. (|6.2I) exists characterized by the following hydrodynamic profiles 

inai: 

~ dT 

p = const, — — = a = const, — — = const. (6.5) 

OS OS 

Note that a = a/ujo\2 oc a/u. Thus, it is of the order of the Knudsen number associated with the 
shear rate. It is important to bear in mind that, since a = const, the ratio a/u is spatially uniform 
even though neither the shear rate a nor the collision frequency u are. Apart from a, there is another 
Knudsen number, this time associated with the thermal gradient. It can be conveniently defined as 

e = V^^. (6.6) 

This quantity is not constant since dgT = const implies e oc T^^/^. As will be seen, the consistency 
of the profiles (16.51) is possible only if a takes a certain particular value for each coefficient of 
restitution a. In contrast, the reduced thermal gradient ? will remain free and so independent of a 

m- 

From Eqs. (16.51) and (16.61) we get 

7^ — = _ = const. (6.7) 

OUx a 
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This means that when the spatial variable {y or s) is eliminated to express T as a function of 
one gets a linear relationship. For this reason, the class of states defined by Eq. (16.51 ) is referred to 
as the LTu class 117511751 [821 . In terms of the variable s, the temperature profile is 



T{s) = To 1 



^0^ 



^^0 



(6.8) 



where To and ?o are the temperature and Knudsen number at a reference point s = 0. 

In order to get the pressure tensor and the heat flux in the LTu flow, it is convenient to introduce 
the dimensionless velocity distribution function 



Tis) 
P 



2T{s) 



m 



d/2 



/(^,v), 



V — US 



v/2T(s)/m' 



(6.9) 



As a normal solution, all the dependence of / on s must occur through the hydrodynamic fields T 
and Ux. Consequently, 

9/ _ dTdf_ ^ du^dl_ 
ds ds dT ds dux' 
Taking into account Eq. (16.91) , and after some algebra, one gets lITSl 



d 



d 



- 2 + — -c + e— 



dc 



81 



dCr 



2{d + 2) 
-0(c)0(ci)] 



dci j dCT [a- V(c')0(c'i^ 



We consider now the following moments of order k = 2r + ^, 

A^2r|V^(?) = j dc^^4-^^0(c;?), Q<h<L 



(6.11) 



(6.12) 



By definition, A^o|o,o = 1> -^o|i,o = A^oji,! = 0, and A^2|o,o = f • According to Eq. (16.111) , the 
moment equations read 



- ( 2r 



1 - e— ) M2r\i+i,h+a. {2rM2r-2\e+2,h+i + hM2r\e,h-i) = J2r\e,h, (6.13) 



where 



J2r\e,h = J dcc2^4-^^:r[c| 



(6.14) 



are the corresponding coUisional moments. In the case of IMM, J'2r\£,h is given as a bilinear 
combinations of the form A^2r'|£',?i'-^2r"|^",/i" such that 2r' + i' + 2r" + i" = 2r + £. Therefore, 
only moments of order equal to or smaller than 2r + i contribute to J2r\i,h- 

It can be verified that the hierarchy (16.13! ) is consistent with solutions of the form 11751 



M 



2r+e-2 
j=0 



fx 



{2r\l,h) 



Oifj + 



odd, 



(6.15) 
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i.e., the moments A^2r-|£,/i(e) of order 2r + £ > 2 are polynomials in the thermal Knudsen number 
e'of degree 2r + £ — 2 and parity i. 

It is important to remark that, when = 0, the hierarchy ( 16.131 ) reduces to that of the (station- 
ary) USF problem for IMM, i.e., the hierarchy obtained from Eq. (15.21) . In other words, the USF 
moments provide the independent terms of the corresponding LTu moments. On the other hand, 
the hierarchy (16.131) reduces to that of the conventional Fourier flow problem for elastic Maxwell 
particles when a = 1 (which, as will be seen below, implies a = 0) [14911701 . The general problem 
(? 7^ 0, a 7^ 0) is much more difficult since it combines both momentum and energy transport. 
The interesting point is that, although the moment hierarchy (16.131) couples moments of order k to 
moments of a higher order A; + 1, it can be exactly solved via a recursive scheme. This is possible 
because the coefficient ^(^^1^+^''^) with j = 2r + ^ — 1 of the moment M.2r\i+i,h do not contribute 
to Eq. (16.131 ). In the following, we will focus on the moments of second order (pressure tensor) 
and third order (heat flux). 

Since the moments of order k are polynomials in e of degree k — 2, it turns out that the ele- 
ments of the pressure tensor are independent of the reduced thermal gradient Therefore, they 
coincide with those obtained in the steady-state USF, being given by Eqs. (15.101 ). (15.141 ). and (15.151 ). 
Moreover, the reduced shear rate a/u is again given by Eq. (15.131 ). so that 



d+1 



a 



d{l + a) 



V2(c/ + 2)(1 



(6.16) 



This confirms that, as said above, the value of a in the LTu flow is enslaved to the coefficient of 
restitution a. If one defines the rheological functions rj and ^ by Eqs. (15.161) and (15.171) . their 
expressions in the LTu flow are the same as those in the USF [see Eqs. (15.181) and (15.191) 1. 

The third-order moments (absent in the USF) are linear functions of ? that cannot be evaluated 
autonomously. However, they can be obtained from Eq. (16.131) in terms of the independent terms 
of the fourth-order moments iTTSl . The explicit forms for the two third-order moments defining the 
X and y components of the heat flux are 



A^2|i,o(e) 



2d? 



U{2d -2 + 3rfC)C/^i°'^'°^ + 8rfC/io°''''^ + (18rf - 18 + 19rfC)/^o''''°^ 



Q\l2dC (2rf - 2 + 3dC)/^i°'^'^^ + Aio 



,(2|2,1) 



(6.17) 



2dC [Ad{ld - 2 + 9rfC)C/io°'^'°^ + 6(6d - 6 + hdQ^f^'^^ 



+3(9rf + 6 + 17dC)yUo ] - 12rfC(7rf - 2 + 9dC)yUo 
-9(6rf - 6 + 5rfC)/i?'^'^^ 



(0|4,1) 



(6.18) 



where 



X = 36(d -ly -d (76 - 56d - 9rfC ) C 



(6.19) 
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and 



c 



C d + 21 



a 



(0|4,1) 



d 1 + a 

,(2|2,1) 



(6.20) 



The coefficients /io°''''°'', /io°'^'^^ /^o^'^'°^ f^o'^'*''^' ^ ^nd /io^'^'^'' are not but the moments A^o|4,05 -^o|4,2. 
-^2|2,05 -^o|4,i5 and M2\2A, respectively, evaluated in the USF [TTSI . This completes the determi- 
nation of A^2|i,o and A^2|i,i- 

Once the non-zero components of the heat flux Qy and are known, one can define an effective 
thermal conductivity n and a cross coefficient $, respectively, by 



Qy 



dT 



— K- 



^ dT durr 



O 7 MX i~\ 

ay a?/ oy 

In the three-dimensional case the expressions of k and $ are 

36 



(4 



$ = $ 



a)2(829 - 162a - 91a2) (7(a) 
864 



o; 



(6.21) 



(6.22) 



(6.23) 



35(1 + a)(4 - a)3(829 - 162a - 91^2) C{a) ' 

where the functions A(a), B(a), and C{a) are polynomials in a of degrees 26, 26, and 24, re- 
spectively, whose coefficients are given in Table 1 of Ref. fTSl . In Eq. (16.231) , $o = |^o'7o/p is the 
corresponding Burnett coefficient in the elastic limit for d = 3 ETl . 



© 
© 

c 




a 



Figure 9: Plot of the reduced heat flux coefficients k/kq (solid line) and $/$o (dashed line) in the 
LTu flow for d = 3. These quantities diverge at a^ — 0.046. 

While the coefficient k is an extension of the conventional NS thermal conductivity, the coef- 
ficient $ is absent at NS order and thus it can be seen as an extension of a Burnett-order transport 
coefficient. Figure |9] depicts the a-dependence of the reduced coefficients k/kq and $/$o- In- 
terestingly, both reduced coefficients are quite similar for the whole range ac < a < 1, where 
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ttc — 0.046, the relative difference being smaller than 10%. As expected, these coefficients di- 
verge when a — )> as a consequence of the similar divergence of the fourth-order moments in the 
USE An important non-Newtonian effect (not directly observed in Fig. [91), is that, as a decreases, 
the magnitude of the streamwise component grows more rapidly than that of the crosswise com- 
ponent Qy, SO that the former becomes larger than the latter for a < 0.9, what represents a strong 
far-from-equilibrium effect. Finally, comparison between the generalized thermal conductivity k 
and the NS coefficient kns shows that, in contrast to the cases of r] and ?7ns5 both coefficients be- 
have in a qualitatively similar way. In the range 0.475 < a < 1, k > (the relative difference 
being smaller than 20%), while k < kns for a < 0.475. 



7. Small spatial perturbations around the USF 

The LTu flow described in the preceding section can be seen as the USF perturbed by the existence 
of a thermal gradient parallel to the velocity gradient (y axis) under the constraints of uniform 
pressure and heat flux. However, the perturbation is not small in the sense that the thermal gradient 
(as measured by the Knudsen number e) is arbitrarily large. In this section, we will adopt a com- 
plementary approach. On the one hand, the perturbations from USF will be assumed to be small, 
but, on the other hand, they will affect all the hydrodynamic fields. 

Let us assume that the USF is disturbed by small spatial perturbations. The response of the 
system to these perturbations gives rise to additional contributions to the momentum and heat 
fluxes, which can be characterized by generalized transport coefficients ||40l . Since the unperturbed 
system is strongly sheared, these generalized transport coefficients are highly nonlinear functions 
of the shear rate. The goal here is to determine the shear-rate dependence of these coefficients for 
IMM. 

To analyze this problem, one has to start from the Boltzmann equation ( 12.11) with a general time 
and space dependence. First, it is convenient to keep using the relative velocity V = v — uq, where 
uq = aySi is the flow velocity of the undisturbed USF state. On the other hand, in the disturbed 
state the true flow velocity u is in general different from uq, i.e., u = uq + 5u, 5u being a small 
perturbation to Uq. As a consequence, the true peculiar velocity is now W = v — u = V — 5u. In 
the Lagrangian frame moving with velocity uq, the Boltzmann equation (12.11 ) reads 

where the gradient V/ in the last term of the left-hand side must be taken at constant V. 

The goal is to find a normal solution to Eq. (17.11) that slightly deviates from the USF. For this 
reason, let us assume that the spatial gradients of the hydrodynamic fields 

A{r,t) = {n{r,t),T{r,t),6u{r,t)} (7.2) 

are small. Under these conditions, it is appropriate to solve Eq. (17.11 ) by means of a generalization 
of the conventional Chapman-Enskog method Il27l . where the velocity distribution function is 



29 



V. Garzo and A. Santos 



Hydrodynamics of inelastic Maxwell models 



expanded about a local shear flow reference state. This type of Chapman-Enskog-like expansion 
has been considered in the case of elastic gases to get the set of shear-rate dependent transport 
coefficients in a thermostatted shear flow problem [|49l |58l and has also been considered in the 
context of inelastic gases [[39ll40ll4n[6Tl . 

As said in sectionHH the Chapman-Enskog method assumes the existence of a normal solution 
in which all space and time dependence of the distribution function occurs through a functional 
dependence on the fields A{v, t), i.e., 

/ = m\A]. (1.3) 

This functional dependence can be made local by an expansion of the distribution function in 
powers of the hydrodynamic gradients: 

/(V) = /(°)(V|A) + e/«(V|A) + ---, (7.4) 

where, as in Eq. (14.21) . e is a bookkeeping parameter that can be set equal to 1 at the end of the 
calculations. The reference zeroth-order distribution function corresponds to the unsteady USE 
distribution function but taking into account the local dependence of the density and temperature 
and the change V — )■ W = V — 5u(r, t). It is important to note that, as seen in section [STTl in 
the stationary USE the temperature is fixed by the shear rate and the coefficient of restitution [cf. 
Eq. (15.131) 1. Therefore, in order to have T as an independent field, one needs to solve the time- 
dependent USE problem, as discussed in section 15.2.1 As a consequence, the associated solution 
f^^\ in dimensionless form, is a function of a and a* = a/u^ separately. Apart from this difficulty, 
a new feature of the Chapman-Enskog-like expansion (17.41 ) (in contrast to the conventional one) is 
that the successive approximations f^'^^ are of order k in the gradients of n, T, and Su, but retain 
all the orders in the reduced shear rate a* ||40| . 

The expansion (17.41) yields the corresponding expansions for the fluxes: 

p = p(o) + ep(i) + . . . , q = eq(i) + . . . , (7.5) 

where Pj^p = pP*j{a,a*) is the pressure tensor in the unsteady USE and we have taken into 
account that q*^*^) = 0. A careful application of the Chapman-Enskog-like expansion to first order 
gives the following forms for the generalized constitutive equations [|40l : 

Pip = -Vijki^e^Uk, (7.6) 

gf ^ = -KijVjT - HijVjU. {1.1) 

In general, the set of generalized transport coefficients rjijM, f^ij, and fiy are nonlinear functions of 
the coefficient of restitution a and the reduced shear rate a*. The anisotropy induced in the system 
by the presence of shear flow gives rise to new transport coefficients, reflecting broken symmetry. 
The momentum flux is expressed in terms of a viscosity tensor r]ijke{a*,a) of rank 4 which is 
symmetric and traceless in ij due to the properties of the pressure tensor Pip . The heat flux is 
expressed in terms of a thermal conductivity tensor Kij{a*,a) and a new tensor Hij(a*,a). Of 
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course, for a* = and a = 1, the usual NS constitutive equations for ordinary gases are recovered 
and the transport coefficients become 

Vijke Vo (^^ik^ji + ^jk^ii - "^^ij^ki^ , i^ij t^oSij, fiij 0. (7.8) 

The elements of the tensor rjij^^ obey a set of coupled linear first-order differential equations 
in terms of a*, P*p and da*Pij. Those differential equations become algebraic equations when one 
specializes to the steady-state condition (|5.6I) . In that case one gets 

a* {SixVjyke + ^jxViyke - ^kyVijxe) + ^o\2Vijke = hia*da*P*j + 6^kP■^ + 5jkP*i 

(7.9) 

where r/*,-;,^ = rjijkt/rjQ. In Eq. (17.91) . a*, P*^, and 5a* P-* are functions of a given by Eqs. (15.111) . 
(I5.13I) - (I5.15I) . (15.271) . and (15.281) . As a simple test, note that in the elastic limit (a* = 0, z/*|2 = 1, 
Pij = Eq. (17.91) becomes Eq. (17.81) . Also, it must be noted that, on physical grounds, the 
elements of the form rjij^y are directly related to the unperturbed transport coefficients rj and \1/ 
defined by Eqs. (15.161) and (15.171) . The rationale is that the particular case of a perturbation in the 
velocity field of the form 5n = {5a)ySi is totally equivalent to an unperturbed USE state with a 
shear rate a + Sa. As a consequence, one has 

P*y{a* + 6a*) = -r]*{a* + 6a*){a* + 6a*) 

= -r]*{a*)a* -6a*{l + a*dl)r]*{a*) + --- , (7.10) 



P^yia* + 6a*) = P:,{a* + 6a*) = 1 - ^^*(a* + 6a*) {a* + 6a*)^ 

1 0^* / n* \ 

= 1 - -**(a*)a*' - ( 1 + y^a J ^*(«*) + ■ ■ ■ , 

(7.11) 

P;,(a* + 6a*) = d-{d- l)P*yy{a* + 6a*). (7.12) 

This implies that 

r]ly,y{a*) = {l + a*d:)r^*{a*), (7.13) 
2fl* / a* \ 

Vlyxyia*) = v:..y{an = ^ + Y^lJ ^*(«*)' (^.14) 

vL.y{an = -id-l)v;yxyia*). (7.15) 

Here, rj* = 77/770, 6a* = 6a/ uq, and = ^/{'ql/p). Equations (I7.13I) - (I7.15I) hold both for the 
steady and unsteady USE. It can be checked that they are consistent with Eq. (17.91 ) in the steady 
state, in which case da*ri* and da*'^* are obtained from Eqs. ( 15.271) and ( 15.281) . The elements of the 
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steady-state shear viscosity tensor ri*jf^^ with k ^ x and 
of algebraic equations (I7.9I ). 



7^ y must be obtained by solving the set 



It turns out that there are two classes of terms. Class I is made of those coefficients r]*jf^^ with 
{k,i) = (xx), (xy), (yx), (yy), and (zz). The complementary class II includes the coefficients 
with {k, £) = (xz), (yz), (zx), and (zy). Of course, class II (as well as the elements r]*j^^ of class 
I) are meaningless if d = 2. The coefficients of the form 77*^^^^ and 77*^^^ vanish in class I, while 
those of the form 77*^^^, Vxyke^ ^^'^ '^yyki vanish in class n. The remaining elements in class II are 
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(7.16) 
(7.17) 

(7.18) 
(7.19) 



Some of the above results might have been anticipated from simple arguments, as shown on p. 138 
ofRef. flU. 

The expressions for the non-zero elements of class I include the derivatives da*r]* and da*"^*. 
Those expressions are much more involved than those of class II and so they will not be explicitly 
given here, except for the cases of Eqs. (I7.13I) - (I7.15I) . As Eq. (17.151) shows, the combination 

Vxxki + (d — l)r]*yi^g vanishes for {k, i) = (xy). It also does for {k, £) = (xx), while for the other 
cases of class I one simply has 



V: 



xxyx 



+ {d- l)r], 



yyyx 



2{d 



p* 

yy 



V*XXyy +{d- l)r]*yyyy = 2 {(1 - 2) , 

Vxxzz ~\~ {d — l)riyy^^ = —2rj . 
Figure [10] shows the steady-state values of two elements of class I {r]. 



(7.20) 

(7.21) 
(7.22) 



xyxy 

and r]ly^y) and two 



elements of class II iri*^zzx ^^d t]*^^^,) in the three-dimensional case. We recall that the first two 



coefficients measure the deviations of Pxy and Pyy with respect to their unperturbed USF values 
due to a perturbation of the form d5ux/dy. Analogously, the two coefficients 'r]*^zzx ^nd "ql^zx 
measure the presence of non-zero values of Pxz and Pyz, respectively, due to a perturbation of the 
form dbuzjdx. We observe that, at a given value of a, the largest influence occurs on P^z- It is 
also interesting to note that rfr^y^y becomes negative at strong values of dissipation, while rfy^zx is 
always negative. 

The evaluation of the heat flux coefficients k.^ and /ijj is more involved than that of the shear 
viscosity tensor ?7jjfc£. In the general unsteady case, and y^j obey coupled linear first-order 
differential equations where, in addition to a*, P^* , and da*Pij, the fourth-order moments of the 
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Figure 10: Plot of the reduced generalized transport coefficients ri^yxy/vo (solid line), riyy^y/riQ 
(dashed line), 7]xzzx/vo (dash-dotted line), and 7]yzzx/vo (dotted line) in the steady-state USF for 
d = 3. 



USF and their first derivatives with respect to a* are also involved [!40l. For steady-state conditions. 



the set of equations becomes algebraic. For further analysis, let us rewrite Eq. (17.71) as 

^ = -KijVjT - JlijVjU + %jVa*, (7.23) 

so that 

In Eq. (|7.23l) we have disentangled the contributions to the heat flux directly associated with 
the temperature and density gradients from those due to the local spatial dependence of a* cx 
^-12^-1/2 Whereas the coefficients Kjj and /ly are given in terms of the second- and fourth-order 
moments of USF, but not of their derivatives with respect to a*, the coefficients are linear func- 
tions of those derivatives. The equations for Kij and /ijj are not reproduced here and the interested 
reader is referred to Ref. [|40l . 

It is illuminating to connect the elements Kij and Jlij in the steady state with the LTu transport 
coefficients k, and $ defined by Eq. (16.211) . The key point is the realization that the LTu (see Sec.l6l 
can be interpreted as a special perturbation of the USF such that (a) the only non-zero temperature 
and density gradients are dyT and dyU, (b) those gradients are not independent but are related by 
the constant-pressure condition dyU = —{n/T)dyT, and (c) the reduced shear rate a* is constant. 
Although in the LTu the strength of the "perturbation" dyT is arbitrary, we have seen that the heat 
flux is linear in the thermal gradient, so that the effective coefficients defined by Eq. (16.211) must 
be related to those defined by (|7.23l) as 

Tl Tl 
K ^yy 7j}l^yyi ^Cl ^xy ~^l^xy C^-^S) 
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We have checked that the above relations are indeed satisfied, what shows the consistency of our 
results. 



8. Concluding remarks 

Exact solutions in nonequilibrium statistical mechanics are scarce. In order to overcome this limi- 
tation, two possible alternatives can be envisaged from a theoretical point of view. On the one hand, 
one can consider a realistic and detailed description but make use of approximate (and sometimes 
uncontrolled) techniques to get quantitative predictions. On the other hand, one can introduce an 
idealized mathematical model (which otherwise captures the most relevant physical properties of 
the underlying system) and solve it by analytical and exact methods. Here we have adopted the 
second strategy. 

In this review we have offered an overview of some recent exact results obtained in the context 
of the Boltzmann equation for a granular gas modeled as inelastic Maxwell particles. Although 
most of the results reviewed in this paper have been reported previously, some other ones are 
original and presented here for the first time. 

We have focused on the hydrodynamic properties of the system, where here the term 'hydro- 
dynamics' has been employed in a wide sense encompassing both Newtonian and non-Newtonian 
behavior. More specifically, the Navier-Stokes (NS) transport coefficients ?7ns, ^^ns, /^ns, and Dns 
have been obtained from the Chapman-Enskog method in Sec. |4] As a common feature, it is ob- 
served that the coUisional inelasticity produces an increase of all the NS transport coefficients, 
especially those related to the heat flux. In fact, the latter coefficients diverge (in two and three 
dimensions) for sufficiently small values of the coefficient of restitution a. This is a consequence 
of the algebraic high- velocity tail of the distribution function in the homogeneous cooling state. 

As an example of non-Newtonian hydrodynamics, the study of the paradigmatic uniform shear 
flow (USF) has been addressed in Sec. [5] The analysis includes both the steady and (transient) 
unsteady states. While the former has been extensively studied in the literature, the hydrodynamic 
transient toward the steady state has received much less attention. We have studied the rheological 
properties (generalized shear viscosity t] and viscometric function ^!), thus assessing the influence 
of inelasticity on momentum transport. Additionally, the a-dependence of the fourth-order velocity 
cumulant 02 in the steady state has been shown. Again, an algebraic high-velocity tail of the USF 
distribution function gives rise to a divergence of 02 for quite small values of a (a < 0.046 for 
rf = 3). 

Next, a more complex state has been analyzed in Sec. |6]] This state (referred to as LTu) is 
actually a class of planar Couette flows characterized by a uniform heat flux, as a consequence of 
an exact balance between viscous heating and coUisional cooling contributions. In contrast to the 
USF (where only momentum flux is present) and to the Fourier flow for an ordinary gas (where 
only heat flux is present), both momentum and heat fluxes coexist in this class of states. It turns 
out that the rheological properties coincide with those of the USF. The most interesting result is 
that the heat flux is exactly proportional to the thermal gradient with coefficients k and $ that are 
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highly nonlinear functions of a. 

Finally, the problem of arbitrary (but small) spatial perturbations to the USF has been con- 
sidered in Sec. IT] Taking the USF as a reference state and carrying out a Chapman-Enskog-like 
expansion about it, a generalized shear viscosity tensor T]ijke and generalized heat-flux tensors Kij 
and are determined. The a-dependence of 'qxyxy, Vyyxy, Vxzzx, and rjy^zx has been explicitly 
given. 

It is worth emphasizing that all the results reviewed in this paper are exact in the context of 
the Boltzmann equation for IMM, regardless of the degree of dissipation. Moreover, all the results 
are explicit, with the exception of those displayed in Fig.[8l which are obtained from a numerical 
solution of the set of coupled differential equations (15.221) and (15.231) . This contrasts with the case 
of inelastic hard spheres, which requires the use of approximations and/or numerical methods. The 
price to be paid is that, in general, the quantitative predictions obtained from IMM significantly 
enhance the influence of dissipation on the dynamical properties of a granular gas. However, 
the IMM is very useful to unveil in a clean way the role played by inelasticity in granular flows, 
especially in situations, such as highly non-Newtonian states, where simple intuition is not enough. 
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